Absence of Critical Thickness in an Ultrathin Improper Ferroelectric Film 
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We study the ferroelectric stability and surface structural properties of an oxygen-terminated 
hexagonal YMnOs ultra-thin film using density functional theory. Under an open circuit boundary 
condition, the ferroelectric state with the spontaneous polarization normal to the (0001) surface, is 
found to be met ast able in a single domain state despite the presence of a depolarizing field. We 
establish a connection between the result and the role of improper ferroelectric transition. Our 
results imply that improper ferroelectric ultrathin films can have rather unique properties that are 
distinctive from those of very thin films of ordinary ferroelectrics. 



The integration of ferroelectric oxide materials into ex- 
isting microelectronic device architectures is currently of 
much interest as they hold promise for a wide range 
of potentially new applications [ll, This has gen- 
erated an enormous effort to understand the proper- 
ties of ferroelectric perovskite ultrathin films, such as 
BaTiOs, where it has been shown that electrical bound- 
ary conditions 0,05 HI, 0], surface and interface prop- 
erties [tI, H, [l, [l^, as well as epitaxial strain [III, 
all play an important role. Although there appears to 
be no fundamental size limit below which ferroelectricity 
disappears, ultrathin ferroelectric films with the polar- 
ization normal to the surface remain a challenge. This 
is because the depolarization field arising from the ac- 
cumulated charges at the surfaces, if not screened, can 
strongly suppress the instability towards a single-domain 
ferroelectric state. This is indeed what happens in a ma- 
terial displaying a proper ferroelectric transition (note, 
henceforth we refer to such materials as proper ferro- 
electrics), namely one where the primary order parameter 
is the electrical polarization, e.g., BaTiOa. The depolar- 
izing field contributes to the free energy a positive term 
quadric in the polarization thereby re-normalizing the 
soft-mode energy and stabilizing the paraelectric phase. 
Even in a system with metallic electrodes, which can 
provide the necessary screening in most cases, when the 
thickness of the ferroelectric film becomes comparable to 
the effective screening length of the metal, the screening 
is incomplete resulting in a reduced polarization [sl. ItI. [i3| . 
an increased coercive field [14], and in some case sup- 
press the tendency towards (single-domain) ferroelectric- 
ity completely. 

It would clearly be advantageous, if not at least funda- 
mentally interesting, to consider the surface properties of 
ultrathin films of materials in which ferroelectricity did 
not originate from a polar instability, but rather from 
an improper ferroelectric transition - i.e., one where the 
spontaneous polarization does not drive the transition 
but instead is a slave to some other primary order param- 
eter. As pointed out decades ago by Levanyuk and San- 
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FIG. 1: Side view of YMnOs (0001) thin film slab two unit- 
cells thick symmetrically terminated with apical oxygens. The 
three Y (large blue circles) planes are indicated. Small red cir- 
cles and polyhedra represent oxygens and MnOs cages. The 
polarization of the unrelaxed FE structure pointing from sur- 
face II to surface I is indicated. 



nikov in a series of papers [lH, [l^ , an instability towards 
a single-domain ferroelectric state in a improper ferro- 
electric is still possible even if the depolarization field 
remains unscreened, e.g., under open circuit boundary 
conditions. 

In this letter we study for the first time finite size ef- 
fects in an improper ferroelectric from first principles. 
We use hexagonal YMnOa [l7|, E, E, 0, which was 
recently predicted [2l| and subsequently experimentally 
confirmed |22|, [13 to display an improper ferroelectric 
transition as an example. We investigate the intrinsic 
properties of isolated thin films of YMnOa with (0001) 
orientation and a clean oxygen terminated surface in the 
absence of an external field or electrode and determine 
the influence of the surface on the ferroelectric state. De- 
spite the growth of YMnOs thin films on a wide range 



of substrates Q , the structural properties of the surface 
have not been extensively studied previously. We show 
that contrary to the case of a thin film proper ferroelec- 
tric, e.g., BaTiOa or PbTiOs 0,11,0, the polarization 
persists in a single domain YMnOs film to approximately 
two unit cells thick (the smallest size in which a bulk re- 
gion can be defined), i.e., the relaxation of the YMnOa 
slab does not result in a paralectric structure. 

We carry out the first-principles study of YMnOs thin 
films using density functional theory as implemented in 
the VASP code [25,] and projected augmented wave pseu- 
dopotentials. The Hubbard correction to the local den- 
sity approximation (LDA+U) is applied to the Mn 3d 
orbitals to account for the correlation effect. We use 
U = 6.0 eV and J = 0.9 eV which was previously 
shown 
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|26| to yield a band gap and lattice param- 
eters close to the experiments. We start with a slab com- 
posed of two unit cells, stacked along the hexagonal axis, 
of the bulk ferroelectric YMnOs structure, terminating 
symmetrically on the apical oxygens with a vacuum layer 
of ~ 10 A added as shown in Fig. [H This is the mini- 
mum thickness for a YMnOs film for which a "bulk" re- 
gion can be defined. We fully relax the atomic positions 
of the slab while keeping the in-plane lattice constants 
fixed at the theoretical bulk ferroelectric equilibrium val- 
ues (i.e., a = 6.09 A and c = 11.36 A). The planewave 
cutoff energy is 500 eV and a Monkhorst-Pack 4x4x1 k- 
point mesh is used for the Brillouin zone integration. We 
have applied collinear approximation to describe the anti- 
ferromagnetic spin ordering while neglecting spin-orbital 
coupling and magnetic anisotropy, which was previously 
found to be small compared to the structural energet- 
ics [nj. 

Previous theoretical work has shown that the chemical 
bonding in rare-earth based hexagonal manganites such 
as YMnOs plays a very minor role in the origin of fer- 
roelectricity [20], unlike that of conventional perovskite 
ferroelectrics where covalency effects between the cations 
and oxygens are crucial [l^ • The difference in the nature 
of bonding from that of perovskites also affects the geo- 
metric relaxation and surface band structure as we will 
show next. 

The geometric relaxation normal to the surface is illus- 
trated in Fig. [21 where I and II denote the same surfaces 
as in Fig. [H At both surfaces, the surface apical (Oa) 
oxygens sitting atop the MnOs cages move substantially 
inward relative to their unrelaxed positions, toward the 
bulk region, with a displacement ranging between 0.25 A 
and 0.33 A. The displacements are larger than those 
reported for, e.g., BaTiOs surface [28], which measure 
of - 0.1 A. This is unsurprising as there are no metal 
cations that bond with the apical oxygens in the com- 
mon surface plane in YMnOs. Instead, the surface oxy- 
gens bond strongly with the Mn atoms in the next atomic 
plane. The surface relaxation should thus be more sig- 
nificant than that found in perovskites. In the next layer 
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FIG. 2: Geometric relaxation of the bottom and top four 
surface atomic planes in the YMnOs slab along the surface 
normal direction. Circles and diamonds depict the unrelaxed 
and relaxed coordinates, respectively, for the surface MnOs 
and Y planes. 



down, the Mn atoms move towards the surface on aver- 
age by 0.03 A and 0.02 A, while the oxygens (Op) move 
towards the bulk region by 0.08 A and 0.17 A for surfaces 
I and II, respectively. The apical oxygens (O^) and Y 
layer below the surface plane have negligible relaxation 
(< 0.01 A). Overall, the interplane distance (calculated 
from position averaged for each layer) at surface I be- 
tween the first two surface planes di2 reduces by 0.23 A 
and between the second and third planes d23 has reduced 
by 0.02A. Similarly, a reduction of di2 by 0.21 A and d23 
by 0.08 A is found at surface 11. 
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FIG. 3: Top view of the in-plane positions of the surface apical 
oxygens (large red circles) relative to the MnOs polyhedra for 
surfaces I and II. The blue medium circles depict Y atoms 
near the surface. The arrows depict the in-plane relaxation 
direction relative to the initial positions. 

Surface I and II differ significantly with respect to their 
in-plane relaxation as shown in Fig. [3l Three apical oxy- 
gens move outward along the in-plane unit-cell axis by 
about 0.38, 0.4, and 0.47 A at surface I, displaying a ten- 
dency to compensate the shortened Mh-Oa bond along 
the c axis. At surface II, two of the three apical oxy- 
gens (in the unit cell) display a tendency to dimerize, 
yielding a 0-0 bond of 2.23 A, 3.14 A, and 3.59 A (the 
length is 3.22 A, 3.67 A, and 3.67 A in the bulk). The 
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FIG. 4: The projected density of states for each MnOs plane 
along the YMnOs thin film c axis as in Fig [1] The top and 
bottom panels correspond to surface II and I. The polarization 
points from the top to the bottom surface. The Fermi level is 
marked at E = 0. 



non-equivalence between the two surfaces evident from 
the in-plane relaxation arises from the presence of the 
ferroelectric polarization as we show below. 

Fig. [4] shows the projected density of states (PDOS) 
for each MnOs plane along the c axis of the film. Like 
the geometric relaxation, the electronic structure of the 
interior layers rapidly converges to that of bulk YMnOa 
where a band gap about 1.4 eV has been reported [29]. 
The top of the valence bands on the surfaces shows more 
changes than the conduction bands with respect to the 
bulk, as the valence band is partially composed of O 2p 
of apical oxygens lying in the surface plane. There are 
surface states that lie deep near the center of the gap 
which pin the Fermi level in the gap. These states are 
primarily of the O 2p character admixed with the Mn 
3d and are localized at the surface plane of the oxygen 
terminated surface. This is somewhat different from a 
perovskite such as BaTiOa, where surface states were 
found to intrude only into the lower part of the gap as a 
result of the strong hybridization between Ti and O [28] . 
Quantitatively, the band structure of the YMnOa film 
relaxed in the FE structure looks similar to that obtained 
for a paraelectric structure (not shown), however, one can 
notice a slight asymmetry in the relaxed FE structure 
between the two surface states, namely, a slight excess 
of electrons at surface I and excess of holes at surface 
II. As we will show the asymmetry can be traced to the 
polarization across the film normal to its surface. 

By determining the potential drop across the fully re- 
laxed YMnOa film in our calculation, we evaluate the 
depolarizing field Ed of 0.025 eV/A across the film [SO^. 
This electric field in the YMnOs slab is approximately 



an order of magnitude smaller than that of an isolated 
PbTiOa slab [7] (the latter can be calculated by arti- 
ficially imposing a constraint on the ions so that the 
structure does not relax back to a paraelectric state) 
and is of course consistent with the fact that the fer- 
roelectric Berry phase polarization [32] in bulk YMnOs 
of 6.5 fiC/cm is approximately one order of magnitude 
smaller than the bulk polarization in PbTiOs. 

The magnitude of the field in a YMnOs film appears 
to indicate that not only does a spontaneous polarization 
exist in an isolated slab of YMnOs but that the value is 
comparable to that in the bulk. Next we calculate, albeit 
approximately, the polarization in the YMnOs slab via 
ionic displacements from a paraelectric reference struc- 
ture to the relaxed ferroelectric structure, u, and Born 
effective charges, Z, where P ^ Z • u [33|. The result- 
ing polarization equals to ~ 6.1 /iC/cm^ which is only 
slightly smaller than the bulk polarization. 

So why does an uncompensated single domain slab of 
BaTiOs (or PbTiOs for that matter) remains paraelectric 
yet a slab of YMnOs under similar conditions becomes 
ferroelectric? The simplest way to understand this dif- 
ference is to investigate a phenomenological free energy. 
In materials that display a proper ferroelectric transi- 
tion such as BaTiOs, the free energy can be expanded in 
terms of the polarization: 

J='{P) = l/2aP^ + l/ApP^ + 0{P^) 

where the ferroelectric transition occurs when a =0 (for 
purposes of illustration we are only discussing second 
order phase transitions). In the case of an uncompen- 
sated film a depolarizing field, Ed oc —P, adds a positive 
term to the free energy, J^d ^ renormalizing a to 
a', where due to the comparatively large depolarization 
energy, a' > at all temperatures, suppressing the fer- 
roelectric transition into a single domain state. The po- 
larization charges on the surface must be compensated 
by free charges provided by, e.g., external electrodes, to 
preserve the polarization. 

In contrast, YMnOs was shown to display an improper 
ferroelectric transition where the polarization arises due 
to a non-trivial coupling to a zone-boundary lattice in- 
stability, the i^-mode. A simplified free energy is given 
by: 

^(P, K) = ao2P^ + a2oK^ + PaoK^ + PsiK^P + ... 

where K corresponds to the primary order parameter as- 
sociated with the zone boundary Ks phonon mode (ac- 
counting for the transition ~1200K) and P is the polar- 
ization [21]. Here the quadratic coefficient of the polar- 
ization, Qfo2, does not soften to zero at the ferroelectric 
transition [3J, i.e., ao2 > at any temperature, and 
the spontaneous polarization arises because of the /^si 
coupling, where for small Pr^ K^. Physically this 
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coupling acts like a field, where Ks "pushes" the equi- 
librium value of the single minimum potential well of 
the polar mode to nonzero value. The energy associated 
with this field, as previously calculated in bulk YMnOa 
(- 5 X 10-"^ eV/A , or 0.06 eV per formula unit), turns 
out roughly 6 times larger than the depolarizing energy. 
For an uncompensated film in this situation, the effect 
of the depolarization field is to stiffen the already " hard" 
quadratic coefficient ao2, which can effectively lower the 
equilibrium value of the polarization (see Ref. |2l|) but 
in principle does not suppress the instability to a single 
domain ferroelectric state. 

We have studied the ferroelectricity and surface relax- 
ation of oxygen terminated YMnOs thin films. Using ul- 
trathin films of YMnOs as an example, we show that in 
improper ferroelectrics, the instability to a single domain 
ferroelectric state does not vanish even when the depo- 
larization field remains unscreened, i.e., in the absence 
of electrodes or external field. Our results are applicable 
for all materials that display an improper ferroelectric 
transition, whether the polarization is induced by struc- 
tural changes, like the case we considered here, or by 
for example a spin ordering. This has potential exciting 
applications given that a wide class of perovskite nanos- 
tructures have been recently found to display improper 
ferroelectricity [IH. 
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